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Summary 


Steady-state, axisymmetric motions of a Boussinesq fluid 
contained in rotating spherical annulus are considered. The 
motions are driven by latitudinally-varying temperature gradients 
at the shells. Linearized formulations for a harrow gap are 
derived and the flow field is divided into the Ekman layers and 
the geostrophic interior. The Ekman layer flows are consistent 
with the known results for cylindrical geometries . Within the 
framework of rather restrictive assumptions, the interior flows 
are solved by a series of associated Legendre polynomials. These 
solutions show qualitative features valid at mid-latitudes . 
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1 . Introduction 

The dynamics of rotating and stratified fluids in a container 
is generally influenced by the presence of Coriolis forces acting 
on the fluid motion in planes perpendicular to the rotation axis. 
Laboratory models have been developed to study the flows to 
simulate the geophysical fluid systems. Historically, cylindrical 
geometries have been widely used because of its maximal simplicity 
and much qualitative information has been made available from 
such model studies . 

One of the most important ingredients in the modeling efforts 
of realistic geophysical fluid systems is the spherically radial 
gravity. Recently, there is work on developing a spherical 
laboratory experiment for Spacelab flights, the Atmospheric 
General Circulation Experiment (AGCE) [see Ref. 1] . Significant 
features of AGCE are that e radial body force , simulating the 
spherical gravity, is imposed via an electric field across 
rotating spherical shells , The fluid contained in the spherical 
gap is stably stratified in the radial direction, and the motion 
is driven by the latitudinal . temperature gradients on the 
boundaries . 

In an effort to obtain qualitative characteristics for such 
experimental configurations, we study the steady, axisymmetric 
motions of a thin layer of fluid betweer two rotating spherical 
shells. As was pointed out in Ref. 2, the axisymmetric flows 
would provide the needed basic-states for stability analysis. 


4 


4 


We shall follow the analytical techniques of Ref, 3, and attempt 
to give an analytical solution within the framework of the 
specified assumptions, The main thrust is in the examination of 
the geostrophic interior flows at mid-latitudes . 

2 , The formulation 

Consider a Boussinesq fluid occupying a spherical annulus 
of thickness D . The radii of the spherical shells , which are 
rotating at angular velocity ft about the polar axis, are 
and R + d , where s s d/r « 1 is assumed, The coordinate 

0 measures latitude, r is the radial distance from the 
sphere's center. The eastward, northward, and radially outward 
velocity components are u, v, w, respectively. 

We introduce the nondimensionalization anticipating that in 
the main body of fluid the flows are in geostrophic and hydrostatic 
balance : 


r=R+6RZ ' , 

(U,V) = uqu'ce.z’), v'(0,z’)], w = <$u w'O.z'), 

'£ = To + AT Z' + (2ftRU/gaD) T' (0 ,Z f ) , 
p = po [1 - a (T - To) ] , 

2 1 

P = Po - p gDZ' + aAT p gD ^ Z r 2 + 2« RU p P'(0 s Z f ), 


for which the usual notation is employed (Ref. 2). 

The equations of motion for linearized, steady-state, 
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axisymmetric motion on a rotating spherical coordinate are (the 
primes have been removed from the nondimensional quantities) i 
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The relevant nondimensional parameters are 

e h = Vh /sir2, e v = Uv/k h 15 = d/r , 
s> = « g AI v D/R 2 R 2 , a H = V H/Kn , a v - V v/K v 

where v v , K v ; v H , Kg are the mixing coefficients for 

momentum and heat in the vertical and horizontal directions, 
respectively. These are introduced to develop a model of greater 
geophysical relevance. 

The boundary conditions are 
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T - t t on Z - 1. (8) 

T » T l on Z « 0. (9) 

(6) and («7) are to be used when there is a mechanical forcing 
of the flow due to a differential rotation of the shells . For 
a purely thermally-driven flow, u L « u T » 0. 

Now, the strategy is to. divide the flow field for mid-latitudes 
into the viscously-controlled Ekman layers on the shells and the 
geostrophic interior . Ref . 2 scales the variables as 


u » 

UjCe.z) + 

u T (e, ; x ) + u. L Ce,? 2 ) . 

(10) 

V = 

E^ce.z) 

+ v T (e,c x ) + v l (6-,5 2 ) . 

(11) 

w = 

w I (6,z) 

+ E v V 8 ’ 5 ! 5 +E v”l < 8 - 

(12) 

p = 

P t (e,z) + 

6 E* 5 ^(0,^) + F l (0,C 2 ). 

(13) 

T = 

T X (S,Z) + 

E v a v S + E V a v S V 8 * «2>- 

(14) 


where the subscript I denotes the interior flow and the tilde 
represents the Ekman layer corrections , c and s are the 
-stretched Z coordinates in the Ekman layers. 

The solutions for the Ekman layer velocities are well-known 
(e.g., Ref.- 4-). The boundary condition for the interior flow is 


TT , n n v ssn0 3 cose[ U I(e,0)- U L] 
w (0,0) - - 2COS0 36 , . 

I sin 0j 


( 15 ) 


and a similar one for (e , 1). 
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conditions for the interior flows are 
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^(0,1) - T t , T x (0,0) » T l . (16) 

3, The geostrophic interior 

At raid-latitudes 0 ~ o (l) , the interior flow equations 

are 
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Following the techniques of Kef. 3 » an integral equation 
involving 3T X / 30 is formed. The solutions of that 

equation may be found in a series of associated Legendre functions 
of order 1 : 
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Prom the thermal wind relation 
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Then, from (20), to o (<S) , we have 

Wj (e,2) - W (O) (0) [(1 + 6) - 2 « Z], 

where w^(0) is given in (22), 
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Now, it is apparent that 
the mid-latitudes e~0(l) so 
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( 29 ) and ( 31) are valid only in 
that sin© ~ 0 (1) . 
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rely heavily on rather severe scaling assumptions* Bine© the 
boundary conditions for the interior flows were based on the 
presence of the Ekman suction, the analysis is applicable to 
rapidly-rotating fluids* The narrow-gap assumption, which is 
essential for deriving the simplified equations of motion, is 
appropriate for realistic atmospheric flows but less accurate 
for laboratory experiments, The interior solutions presented 
lose validity as the equator is approached, The horizontal 
component of the rotation vector vanishes at the equator, and, 
therefore, the Ekman layer fades away, The interior flows 
themselves should satisfy the boundary conditions on the shells, 
and the viscously-controlled horizontal velocities near the equator 
have to be accounted for. To remove the singularity at the 
equator, an equatorial boundary layer is called for. Analytical 
treatments of this equatorial boundary layer appear to be very 
difficult and the attempts so far have been inconclusive. The 
results of the interior flows for e~o (1) show features 
qualitatively consistent with other studies. A full-length paper 
is being prepared to report on the findings. 
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Summary 

Steady-state, axisymmetric motions of a Boussinesq fluid 
contained in rotating spherical annulus are considered, The 
motions are driven by latitudinally-varying temperature gradients 
at the shells , Linearized formulations for a narrow gap are 
derived and the flow field is divided into the Ekman layers and 
the geostrophic interior. The Ekman layer flows are consistent 
with the known results for cylindrical geometries. Within the 
framework of rather restrictive assumptions, the interior flows 
are solved by a series of associated Legendre polynomials, 
solutions show qualitative features valid at mid-latitudes . 
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1 . 


Introduction 


The dynamics of rotating and stratified fluids in a container 
is generally influenced by the presence of Coriolis forces acting 
on the fluid motion in planes perpendicular to the rotation axis, 
Laboratory models have been developed to study the flows to 
simulate the geophysical fluid systems. Historically, cylindrical 
geometries have been widely used because of its maximal simplicity 
and much qualitative information has been made available from 
such model studies , 

One of the most important ingredients in the modeling efforts 
of realistic geophysical fluid systems is the spherically radial 
gravity. Recently, there is work on developing a spherical 
laboratory experiment for Spacelab flights , the Atmospheric 
General Circulation Experiment (AGCE) [see Ref. 1] . Significant 
features of AGCE are that a radial body force , simulating the 
spherical gravity, is imposed via an electric field across 
rotating spherical shells. The fluid contained in the spherical 
gap is stably stratified in the radial direction, and the motion 
is driven by the latitudinal . temperature gradients on the 
boundaries . 

In an effort to obtain qualitative characteristics for such 
experimental configurations, we study the steady, axisymmetric 
motions of a thin layer of fluid between two rotating spherical 
shells. As was pointed out in Ref. 2, the axisymmetric flows 
would provide the needed basic-states for stability analysis. 



We shall follow the analytical techniques of Ref, 3, and attempt 
to give an analytical solution within the framework of the 
specified assumptions . The main thrust is in the examination of 
the geostrophic interior flows at mid-latitudes . 

2 . The formulation 

Consider a Boussinesq fluid occupying a spherical annulus 
of thickness D , The radii of the spherical shells, which are 
rotating at angular velocity n about the polar axis, are 
and R + D , where $ = D/R « 1 is assumed. The coordinate 

6 measures latitude, r is the radial distance from the 
sphere's center. The eastward, northward, and radially outward 
velocity components are U, v, w, respectively. 

We introduce the nondimensionalization anticipating that in 
the main body of fluid the flows are in geostrophic and hydrostatic 
balance : 


r=R+5RZ' , 

(U,v) = U [u' (0 ,z r ) , v r (0 ,z ' ) ] , w = csu W'(e,z'), 

T = To + AT v Z' + (2ftRU/gaD) T'(6,Z r ), 

P = po [1 - a (T - To) ] , 

2 i _ 

P = Po - P o gDZ' + aAT v P Q gD j Z'2 + 2fl RU P o P'CQ,Z T ), 


for which the usual notation is employed (Ref. 2). 

The equations of motion for linearized, steady-state, 
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The relevant nondime ns ional parameters are 
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where v v » K v ? V H > ^ are the mixing coefficients for 

momentum and heat in the vertical and horizontal directions, 
respectively. These are introduced to develop a model of greater 
geophysical relevance , 
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T - t t on z - 1. (8) 

T » T l on Z « 0. (9) 

(fc) and (7) are to be used when there is a mechanical forcing 
of the flow due to a differential rotation of the shells . For 
a purely thermally-driven flow, u L « u T « 0. 

Now, the strategy is to. divide the flow field for mid-latitudes 
into the viscously-controlled Etonian layers on the shells and the 
geostrophic interior. Ref. 2 scales the variables as 
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where the subscript I denotes the interior flow and the tilde 
represents the Ekman layer corrections , c and c are the 
E ^ -stretched z coordinates in the Ekman layers . 

The solutions for the Ekman layer velocities are well-known 
(e.g., Ref 4) . The boundary condition for the interior flow is 
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The correction to the temperature in the EAman layer is 
o(a v S Ep } thus, to OCeJj ) , the temperature boundary 

conditions for the interior flows are 
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3 , The geostrophic interior 


At mid-latitudes 0-0 (l) , the interior flow equations 
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Following the techniques of Ref, 3, an integral eqtiation 
involving 3T X / 36 is formed, The solutions of that 

equation may be found in a series of associated Legendre functions 
of order 1 > 
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Then, from (20), to Q(6) , we have 
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Now, it is apparent that 
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( 29 ) and ( 3l) are valid only in 
that sin0 ~ 0 (1) . 


The above analysis for the interior flows at mid-latitudes 


rely heavily on rather severe scaling assumptions. Since the 
boundary conditions for the interior flows were based on the 
presence of the Ekman suction, the analysis is applicable to 
rapidly-rotating fluids, The narrow-gap assumption, which is 
essential for deriving the simplified equations of motion, is 
appropriate for realistic atmospheric flows but less accurate 
for laboratory experiments. The interior solutions presented 
lose validity as the equator is approached. The horizontal 
component of the rotation vector vanishes at the equator, and, 
therefore, the Ekman layer fades away. The interior flows 
themselves should satisfy the boundary conditions on the shells , 
and the viscously-controlled horizontal velocities near the equator 
have to be accounted for. To remove the singularity at the 
equator, an equatorial boundary layer is called for. Analytical 
treatments of this equatorial boundary layer appear to be very 
difficult and the attempts so far have been inconclusive . The 
results of the interior flows for e~o (i) show features 
qualitatively consistent with other studies. A full-length paper 
is being prepared to report on the findings, 
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